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Abstract —Achieving long battery lives or even self sustain¬ 
ability has been a long standing challenge for designing mobile 
devices. This paper presents a novel solution that seamlessly in¬ 
tegrates two technologies, mobile cloud computing and microwave 
power transfer (MPT), to enable computation in passive low- 
complexity devices such as sensors and wearable computing de¬ 
vices. Specifically, considering a single-user system, a base station 
(BS) either transfers power to or offloads computation from a 
mobile to the cloud; the mobile uses harvested energy to compute 
given data either locally or by offloading. A framework for energy 
efficient computing is proposed that comprises a set of policies 
for controlling CPU cycles for the mode of local computing, 
time division between MPT and offloading for the other mode 
of offloading, and mode selection. Given the CPU-cycle statistics 
information and channel state information (CSI), the policies aim 
at maximizing the probability of successfully computing given 
data, called computing probability, under the energy harvesting 
and deadline constraints. The policy optimization is translated 
into the equivalent problems of minimizing the mobile energy 
consumption for local computing and maximizing the mobile 
energy savings for offloading which are solved using convex 
optimization theory. The structures of the resultant policies are 
characterized in closed form. Furthermore, given non-causal CSI, 
the said analytical framework is further developed to support 
computation load allocation over multiple channel realizations, 
which further increases the computing probability. Last, simu¬ 
lation demonstrates the feasibility of wirelessly powered mobile 
cloud computing and the gain of its optimal control. 

Index Terms —Wireless power transfer, energy harvesting com¬ 
munications, mobile cloud computing, energy efficient computing. 


I. Introduction 

The explosive growth of Internet of Things (IoT) and mobile 
communication is leading to the deployment of tens of bil¬ 
lions of cloud-based mobile sensors and wearable computing 
devices in near future [ lj. Prolonging their battery lives and 
enhancing their computing capabilities are two key design 
challenges. They can be tackled by several promising tech¬ 
nologies; 1) microwave power transfer (MPT) for powering 
the mobiles using microwaves 2) mobile computation 
offloading (MCO) for offloading computation-intensive tasks 
from the mobiles to the cloud 0, and 3) energy efficient 
local computing using the mobile CPU. These technologies 
are seamlessly integrated in the current work to develop 
a novel design framework for realizing wirelessly powered 
mobile cloud computing under the criterion of maximizing 
the probability of successfully computing given data, called 
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computing probability. The framework is feasible since MPT 
has been proven in various experiments for powering small 
devices such as sensors (4j or even small-scale airplanes and 
helicopters 0. Furthermore, sensors and wearable computing 
devices targeted in the framework are expected to be connected 
by the cloud-based IoT in the future j l j, providing a suitable 
platform for realizing MCO. 

A. Prior Work 

MCO has been an active research area in computer sci¬ 
ence 0 where research has focused on designing mobile- 
cloud systems and software architectures |6|. (7j|, virtual 
machine migration design in the cloud |8) and code parti¬ 
tioning techniques in the mobiles [6 j for reducing the energy 
consumption and improving the computing performance of 
mobiles. Nevertheless, implementation of MCO requires data 
transmission and message passing over wireless channels, 
incurring transmission power consumption (9J. The existence 
of such a tradeoff has motivated cross-disciplinary research on 
jointly designing MCO and adaptive transmission algorithms 
to maximize the mobile energy savings fTO) -m A stochastic 
control algorithm was proposed in [T()| for adapting the 
offloaded components of an application to a time-varying wire¬ 
less channel. Furthermore, multiuser computation offloading in 
a multi-cell system was explored in jjllj, where the radio and 
computational resources were jointly allocated for maximizing 
the energy savings under the latency constraints. In p~2) , the 
threshold-based offloading policy was derived for the system 
with intermittent connectivity between the mobile and cloud. 

Energy-efficient mobile (local) computing is also an active 
field where rich and diversified techniques have been designed 
for reducing the mobile energy consumption The 

scheduling of multiple computing tasks was optimized in 
G3 to increase the energy savings. For the same objective, 
dynamic power management was proposed in 0 where com¬ 
ponents of computing tasks were dynamically reconfigured 
and selectively turned off. Another energy-efficient approach 
is to control the CPU-cycle frequencies under the deadline 
constraint by exploiting the fact that lowering the frequencies 
reduces the energy consumption |15|-jT7]|. Recently, the CPU- 
cycle frequencies are jointly controlled with MCO given a 
stochastic wireless channel in |T8). The framework is further 
developed in the current work to include the new feature of 
MPT. This introduces several new design challenges. Among 
others, the algorithmic design of local computing and of¬ 
floading becomes more complex under the energy harvesting 
constraint due to MPT, which prevents energy consumption 
from exceeding the amount of harvested energy at every time 
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instant G3- Another challenge is that MPT and offloading 
time shares the mobile antenna and the time division has to 
be optimized. 

The MPT technology has been developed for point-to- 
point high power transmission in the past decades [Zj. Re¬ 
cently, the technology is being further developed to power 
wireless communications. This has resulted in the emergence 
of an active field called simultaneous wireless information 
and power transfer (SWIPT). Various SWIPT techniques 
have been developed by integrating MPT with communication 
techniques including MIMO transmission ]20], OFDMA [21], 
two-way communication [ 22] and relaying |23|. Furthermore, 
existing wireless networks such as cognitive radio |24) and 
cellular networks |25) have been redesigned to feature MPT. 
Recent advancements on SWIPT are surveyed in | [26) and |27| . 
Most prior work on SWIPT aims at optimizing communica¬ 
tion techniques to maximize the MPT efficiency and system 
throughput. In contrast, the current work focuses on optimizing 
the local computing and offloading under a different design 
criterion of maximum computing probability. 


B. Contributions and Organization 

Consider a single-user system comprising one multi-antenna 
base station (BS) using transmit/receive beamforming for 
transferring power to a single-antenna mobile or relaying of¬ 
floaded data from the mobile to the cloud. To compute a fixed 
amount of data, the mobile operates in one of the two available 
modes: local computing and offloading. In the mode of local 
computing, MPT occurs simultaneously as computing based 
on the controllable CPU-cycle frequencies. Nevertheless, in 
the mode of offloading, the given computation duration is 
adaptively partitioned for separate MPT and offloading since 
they share the mobile antenna. Assume that the mobile has the 
knowledge of statistics information of CPU cycles and channel 
state information (CSI). The individual modes as well as mode 
selection are optimized for maximizing the computing prob¬ 
ability under the energy harvesting and deadline constraints. 
For tractability, the metric is transformed into equivalent ones, 
namely average mobile energy consumption and mobile energy 
savings, for the modes of local computing and offloading, 
respectively. Compared with ||T8|, the current work integrates 
MPT with the mobile cloud computing, which introduces 
new theoretical challenges. In particular, the energy harvesting 
constraint arising from MPT makes the optimization problem 
for local computing non-convex. To tackle the challenge, the 
convex relaxation technique is applied without compromising 
the optimality of the solution. It is shown in the sequel that 
the local computing policy of [18] is a special case of the 
current work where the transferred power is sufficiently high. 
Furthermore, the case of dynamic channel for mobile cloud 
computing is explored. Approximation methods are used for 
deriving the simple and close-to-optimal policies. 

The contributions of the current work are summarized as 
follows. 

• Optimal local computing: First, consider a static channel. 

For the mode of local computing, CPU-cycle frequencies 


are optimized for minimizing the average mobile energy 
consumption under the energy harvesting and deadline 
constraints. The corresponding optimization problem is 
non-convex but solved by the convex relaxation that is 
proved to maintain the optimality. The resultant policy 
is shown to depend on two derived thresholds on the 
transferred power. The first determines the feasibility 
of successful computing. Given feasible computing, the 
optimal frequency of the n-th CPU cycle is shown to be 
proportional to p k 3 with pk being the probability of its 
occurrence and independent of the transferred power if it 
is above the second threshold; otherwise, the frequency 
is proportional to (p^ + A) - s with A being a constant 
determined by the transferred power. 

• Optimal computation offloading : Next, for the mode of 
offloading, the partition of the computation duration for 
separate MPT and offloading is optimized to maximize 
the energy savings. As a result, a threshold is derived on 
the product of the BS transmission power and squared 
channel power gain, above which offloading is feasible. 
Given feasibility, the optimal offloading duration is shown 
to be proportional to the input data size and inversely pro¬ 
portional to the channel bandwidth. In other words, small 
data size and large bandwidth reduces time allocated for 
offloading and increases time for MPT and vice versa. 

• Mobile mode selection: The above results are combined 
to select the mobile mode for maximizing the computing 
probability. Given feasible computing in both modes, the 
one yielding the larger energy savings is preferred and 
the selection criterion is derived in terms of thresholds 
on the BS transmission power as well as the deadline for 
computing. 

• Optimal data allocation for a dynamic channel: Last, the 
above results are extended to the case of a dynamic chan¬ 
nel, modeled as independent and identically distributed 
(i.i.d.) block fading, and non-causal CSI at the mobile 
(acquired from e.g., channel prediction). The problem of 
optimizing an individual mobile mode (local computing 
or offloading) is formulated based on the master-and- 
slave model using the same metric as the fixed-channel 
counterpart. The master problem concerns the optimal 
data allocation for computing in a fixed number of fading 
blocks. The slave problem targets the mode optimization 
in a single fading block for allocated data and is similar 
to the fixed-channel counterpart. By approximating the 
master problems, sub-optimal policies are designed for 
data allocation and shown by simulation to be close-to- 
optimal. The results can be straightforwardly combined 
to enable the mode selection. 

The remainder of this paper is organized as follows. The 
system model is introduced in Section II. Section III presents 
the optimal policies for mobile mode optimization and selec¬ 
tion for the case of static channel. The results are extended 
in Section IV to the case of dynamic channel. Simulation 
results are presented in Section V, followed by the conclusion 
in Section VI. 
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(a) Wirelessly powered mobile cloud computing system 



t’ T time 


(b) Mobile operation modes 

Figure 1: (a) Wirelessly powered mobile cloud computing system and (b) the mobile operation modes. 


II. System Model 

Consider a wirelessly powered mobile cloud computing 
system model shown in Fig. |TJa) comprising one single¬ 
antenna mobile, e.g., mobile sensor and wearable computing 
device with milliwatt power consumption, and one multi¬ 
antenna BS that is a node of a cloud. The BS either transfers 
power wirelessly to or offloads computation from the mobile. 
Consider the computation of a single task that cannot be split 
for partial local computing and partial offloading. In other 
words, the mobile operation mode is either local computing 
or offloading but not both (see Fig. [T] b)). The local com¬ 
puting and MPT can be performed simultaneously while the 
offloading and MPT are non-overlapping in time assuming 
half-duplex transmission. The current work can be generalized 
to the case of computing a multi-task program which allows 
program partition and thus the simultaneous operation of 
these two modes. The multi-antenna BS uses beamforming 
to transfer power or receive the signal. Assuming channel 
reciprocity ©, (21), the effective scalar channel power gain 
is represented as h. Both cases of static and dynamic channel 
are considered. The dynamic channel is modeled as i.i.d. block 
fading where the channel power gain is fixed in each fading 
block and i.i.d. over different blocks. Let P[> denote the BS 
transmission power. It is assumed that both I), and h are 
accurately estimatecQ by the mobile and used for controlling 
the local computing and offloading. During MPT, assuming the 

'Accounting for imperfect channel knowledge, robust optimization tech¬ 
niques can be applied to the current framework for obtaining control policies 
to guarantee a required computing probability. 


mobile has infinite battery storage]^ since the energy harvested 
from the noise is negligible, the energy harvested by the 
mobile in one time unit is v I\h where the constant 0 < v < 1 
represents the energy conversion efficiency. 


A. Local Computing Model 

Adopting the model in fl7) , [18], the required number of 
CPU cycles for computing given data is randomly generated 
as follows. Let L denote the number of bits in input data 
for computation and T represent its deadline. Define the CPU 
cycle information (CCI) as the number of CPU cycles required 
for computing 1-bit of input data. Then the CCI can be 
modeled as a random variable denoted as X and LX gives the 
number of CPU cycles for computing /.-bit input data (1§. 
Define No as a positive integer such that Pr(X > No) < e 
where e is a small real number: 0 < e -C 1. It follows that 
Pr (LX > N) < e where N = LNq. Then given L-bit 
input data, N upper bounds the number of CPU cycles almost 
surely. Define = Pr(LA' > k ) with k = 1,2, ••• ,N 
such that pt, is the probability that the data has not been 
processed completely after k CPU cycles. In other words, pk 
is the probability that the fc-th cycle is executed. Note that the 
sequence pi,P 2 , m " ,pn is monotone decreasing, regardless 
of any probability distribution of the random variable X . 

2 Considering finite battery storage at mobiles will complicate the corre¬ 
sponding structure of the optimal solution. However, it is expected to have 
the threshold-based structure as we derive for the case of infinite storage. 
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Assumption 1. The mobile only has the knowledge of the 
CCI distribution^] 

The energy consumption of local computing is modeled as 
follows. In practice, mobile energy consumption is contributed 
by computation, transmission and fixed circuit power. The 
circuit power is omitted for simplicity but it can be accounted 
for by adding a constant in the problem formulation that, 
however, does not affect the solution method and key results. 
Let E(f) denote the amount of energy for computation in a 
single CPU cycle with the frequency / (or equivalently having 
a CPU-cycle duration of 1//). Following the model in [281, 
p9| , under the assumption of a low CPU voltage, E(f) - 7 f 2 
where 7 is a constant determined by the switched capacitance. 
However, the analysis of current work also applies to a general 
polynomial function E(f) that is monotone increasing and 
convex. Let f \. fa , • • • . /..v denote the CPU-cycle frequencies 
for CPU cycles 1, 2, • • • , TV, respectively. 

Last, for tractability, several assumptions are made for the 
case of dynamic channel. 

Assumption 2. For the case of dynamic channel, the mobile 
has non-causal CSlQ 

Non-causal CSI in Assumption [2] corresponds to the case 
where the mobile predicts the channel variation or has a pre¬ 
determined trajectory. The assumption allows off-line data 
allocation for computing in multiple fading blocks. 

Assumption 3. For the case of dynamic channel, the number 
of channel realizations in the computation duration [0, T] is 
much smaller than the typical number of CPU cycles required 
for mobile local computing in the same duration. 

Assumption [3] states different time scales for the channel 
variation and CPU cycles, allowing energy-efficient control 
policies to be designed using a decomposition approach in the 
sequel. 

Assumption 4. For the case of dynamic channel, the input 
data can be divided continuously into parts that can be 
computed separately. 

The input data of several applications can be divided, such 
as Gzip compression and feature extraction. Assumption [4] 
made for tractability simplifies the data allocation for comput¬ 
ing in separate fading blocks to reduce energy consumption. 
In practice, the optimal solution in the current paper can be 
discretized by rounding. 

B. Computation Offloading Model 

Consider computation offloading where the mobile transmits 
the data to the BS for computing in the cloud and receives the 
result via the BS. Given the mobile transmission power P t , 

3 The results on CPU-cycle optimization in the current work can be straight¬ 
forwardly extended to the case where the mobile has the prior knowledge of 
CCI instead of CCI distribution, resulting in a simpler policy structure. 

Considering causal CSI, if the channel gains of adjacent fading blocks are 
temporally correlated that can be modeled as the channel state Markov chain, 
dynamic programming can be used for deriving the optimal solution, which, 
however, is intractable to analyze and thus is not considered in the current 
work. 


the uplink channel capacity (in bit/s), denoted as C, is given 
by: 


C = Blog 




where B is the channel bandwidth and a 2 is the variance 
of complex white Gaussian channel noise. Moreover, it is 
assumed that the time for computing in the cloud and trans¬ 
mitting the computation result from the BS to the mobile 
is negligible since the cloud has practically infinite compu¬ 
tational resources and the BS can afford high transmission 
power to reduce the downlink transmission delay. Last, the 
computation result is assumed of a small size such that 
demodulating the result data at the mobile results in negligible 
energy consumption compared with that for local computing 
or offloading. 


C. Performance Metrics 

Given the CCI A', let -Empt(() represent the amount of 
cumulative energy harvested by the mobile over duration [0, f] 
and E mob (t, X, {fk} k =i) as the amount of cumulative energy 
consumed by the mobile for local computing or offloading. 
Then the performance metric, computing probability, which 
is the probability of successfully computing given data, is 
denoted as P c , defined as 


Pc= E [I (E M Mt)>E mob (t,X, {fk}k=i), V f € [0,T])] 

(1) 

where the indicator function I(£) gives 1 when the event 
£ occurs and 0 otherwise. Note that the argument of the 
indicator function in ([T} specifies the energy harvesting con¬ 
straint. It can be observed from <0 that maximizing the com¬ 
puting probability by optimizing the CPU-cycle frequencies 
is equivalent to minimizing the mobile energy consumption 
E mo b{t, X, {fk}k=i) f° r ^ £ [0, r]. Nevertheless, this is 
infeasible since X is unknown to the mobile based on Assump¬ 
tion [T] To overcome this difficulty, the current work instead 
focuses on minimizing the expected energy consumption in 
the duration [0, T\, namely E [L/ mo b(T', A', {fk} k =ii] > but still 
under the same constraints, which is given by 

N LX 

E[E moh (T,X,{f k }X =1 )]=Y,JPkti+ E TPfc/k- (2) 

fc=l k=N+1 


Furthermore, for a small e that is close to zero (see Sec¬ 
tion [TLAi, the second component of ([2} is negligible. Then 
the following approximation is used for tractability: 


N 


E [E moh (T, X, {/ fe }f =1 )] « E W*- 

fc=1 


(3) 


Next, consider offloading. By slight abuse of notation, let 
EMPr{t') represent the amount of total energy harvested by 
the mobile over duration [0,/'] and E mob (t') as the amount of 
energy consumed by the mobile for offloading in the duration 
[f, 7’]. Then the computing probability for the current case 
can be defined as 

P C =E[I (E UPT (f) > E mob {t'), Vf'6 [0, T])]. (4) 
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As observed from 0- maximizing the computing prob¬ 
ability is equivalent to maximizing the energy savings 
[-E'MPtV 7 ) — -E-mob^ 7 )]- 

III. Energy Efficient Mobile Cloud Computing 
with a Static Channel 

In this section, given a static channel, the CPU-cycle fre¬ 
quencies and MPT-and-offloading time division are optimized 
for local computing and offloading, respectively. Then the 
results are combined for optimizing the mobile mode selection. 


A. Energy Efficient Local Computing with a Static Channel 

1) Problem Formulation: Based on the discussion in Sec¬ 
tion II-C the problem of optimizing CPU-cycle frequencies 


aims at minimizing the average mobile energy consumption in 
(|3]l under two constraints. The first is the deadline constraint: 
J2k =l 77 — 7 ■ The second is the energy harvesting constraint 
comprising of N sub-constraints given as 


m m 

X^T/fc < vP b hJ2y’ 

k= 1 , ' fc 


m = 1,2,' 


,JV 


(5) 


fc=t 


where the left-hand side of the inequality is the total energy 
consumed by the first m CPU cycles and the right-hand side is 
the total energy harvested by the end of m-th cycle. It follows 
that the optimization problem is formulated as: 


N 


mm 

{/*} 


s.t. 


^ZlVkf\ 


k =1 


ifk < vP b h J2 


k =1 
N 


k =1 


fk 


m = 1, 2, ■ 


•• ,N, 


(PI) 


^TP t ' 

fe=1 

fk > o, 


k = 1,2,- ,N. 


Note that pk is included in the objective function to formulate 
the average energy consumption, while the energy harvesting 
constraint is defined for all possible CPU-cycle realizations 
and thus is without pk- 

2) Solution: It can be observed that the energy harvesting 
sub-constraints in Problem PI are non-convex, resulting in a 
non-convex optimization problem. To address this issue, first, 
define a set of new variables {yk} as yk = jj; for all k. 
Substituting them into Problem PI and furthermore relaxing 
the equality constraint ykfk = 1 to be ykfk > 1 gives 


min 

{fk,Vk} 


s.t. 


N 

J^TPkfk 

k =1 

m m 

Vk, 

k =1 k =1 

N 

^2vk<T, 

k =1 

fk > 0, jr ~yk< o, 


rn= 1,2, ••• ,N, 

(P2) 


k= 1,2, - - - ,N. 


Observed that Problem P2 is a convex optimization problem. 
Nevertheless, the relaxation mentioned earlier has no effect on 
the solution optimality as shown in the following lemma. 

Lemma 1. The solution for Problem P2 also solves PI. 


Proof: See Appendix [A] □ 

Lemma[T]is essential for solving the non-convex Problem PI 
by equating it with the convex Problem P2 that yields the same 
solution but admits powerful convex optimization algorithms. 

To characterize the structures of the optimal CPU-cycle fre¬ 
quencies, several useful properties of the solution for Problem 
P2 are given as follows. 


Lemma 2. The solution for Problem P2, denoted as { yf. ff } , 
satisfies the following. 

1) The deadline constraint is active: 

■spN * _ 'spN J_ _ rp 

Zjfc=l Uk — 2^fc=l /* — ± ■ 

2) The optimal CPU-cycle frequencies can be written as 


p ~ 

vP b h 


)1 

27 

Pk + Em=t 



Vfc (6) 


where the nonnegative variables p, {A m } are the La¬ 
grange multipliers associated with the deadline and en¬ 
ergy harvesting constraints, respectively. 

3) Furthermore, ft </!••■</£. 


Proof: See Appendix [B] □ 

Another important property is stated in the following lemma. 


Lemma 3. The solution for Problem P2 can also be derived by 
solving the following Problem P3 that results from removing 
the first (N — 1) energy harvesting sub-constraints of P2: 

N 


min 

{fk,Vk} 


s.t. 


TPfc/fe 

k =1 

N N 

£ 7/2 < vP b hJ2y k , 

k =1 k—1 

N 

^2 y k < 

*:=1 


(P3) 


fk > 0, jr-yk< 0, k = l,2,---,N. 

Proof: See Appendix [C] □ 

Lemma [3] means the Lagrange multipliers associated with 
the first ( N— 1) energy harvesting sub-constraints are equal to 
zero, i.e., Afc = 0 for k = 1, 2, • • • N — 1. Combining Property 
2) in Lemma [2] and Lemma [3] simplifies the expression for the 
optimal CPU-cycle frequencies as 


Jk 


p — vP b h\ 3 
27 (Pk + A)_ 


Vfc 


(7) 


where A..y is re-denoted as A for ease of notation. To obtain 
the closed-form expressions for {/£}, dehne two positive 
constants a and a' as 


7(V 3 , 

a = —=pr and a 
vl d 


7 

vT 3 




( 8 ) 
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Then the main result of this subsection is stated as follows. 


Theorem 1. (Optimal CPU-cycle frequencies for local com¬ 
puting). The optimal CPU-cycle frequencies {/*, , ■ • • , ff} 

that solve the optimization problem P3 satisfy the following. 

1) If Pbh < a, {/*, / jV '' j/jvl i s an empty set since 
Problem P3 is infeasible. 

2) If a < Pbh < a'. 


fk = 


1 N 


A); 


(Pk + A) *, Vfc (9) 


m —1 

where the positive constant A is the Lagrange multiplier 
with respect to (w.r.t) the energy harvesting constraint of 
P3 and it satisfies 

" N 

J2^ Pk + A) 3 
_fc=l 

3) If Pbh > a! , {/£} are independent of Pbh , 


N 




.*:=i 


vP b hT 3 

7 


( 10 ) 


/ fc * = (^E^Jpfc 4 . Vfc - (id 

\ m=l / 

Proof: See Appendix |D| □ 

The first case ( Pbh < a) corresponds to the scenario where 
the transferred power Pbh is so low that it is impossible to 
complete the local computing within the deadline T. In the 
second case (a < Pbh < a'), the transferred power is not 
large but sufficient for meeting the deadline and consequently 
the optimal CPU-cycle frequencies are functions of Pbh. The 
transferred power is large in the last case (Pbh > a') where 
increasing Pb no longer has any effect on the optimal CPU- 
cycle frequencies and only increases the amount of energy 
savings. In this case, the optimal CPU-cycle frequencies 
depend only on the CCI distribution and the deadline. 


Remark 1. As observed from the definition of a in and 
Theorem [T] there exists a tradeoff between the deadline T 
and BS transmission power Pb. Specifically, meeting a stricter 
deadline requires larger Pb and vice versa. In other words, 
the BS can control Pb to increase the computing probability. 
However, jointly designing the control policies at the BS and 
mobile is challenging and Pb is assumed fixed for simplicity. 


Remark 2 (BS transmission power and the parameter A). 
Recall that A is the Lagrange multiplier for solving Problem 
P3. The optimal CPU-cycle frequencies in Cases 2) and 3) of 
Theorem [I] correspond to A > 0 and A = 0, respectively. The 
case of A = 0 (or equivalently Pbh>a') implies that the energy 
harvesting constraint in Problem P3 is inactive at the solution 
point, reducing the problem and solution to be identical to 
those in fl8) considering local computing without MPT. 


Corollary 1 (Minimum average energy consumption). It 
follows from Theorem |T| that the minimum average energy 
consumption for the local computing, denoted as Ef oc , is given 
as follows. 

1) If a < Pbh < a'. 


F* 

-'HOC rp 2 


E(pfc + A )" 

_k =1 


" N 

_fc= i 


( 12 ) 


and Lj* c is a monotone-decreasing function of Pbh satis¬ 
fying 


1_ 

J^2 



< Kc < 


J^2 


N 

■X>- 

k =1 


2) If P b h > a', 


77'* _ 

Aoc _ 


rp2 



(13) 


that is independent of Pbh. 

Moreover, the corresponding maximum average mobile energy 
savings, denoted as S* oc , is given as S^ oc = vPbhT — E* oc . 

Proof: See Appendix [E] □ 


B. Energy Efficient Offloading with a Static Channel 

This sub-section focuses on computation offloading. The 
time division between MPT and offloading is optimized for 
maximizing the mobile energy savings (see Section II-C| ). 

1) Problem Formulation: The objective function, namely 
the mobile energy savings, is obtained as follows. As shown 
in Fig. [TJb), for current operation mode and t' £ (0 ,T), 
the time interval [0, T\ is divided into two parts: [0, t'} and 
corresponding to MPT and offloading, respectively. 
Let the amount of energy harvested over the interval [0, tj] 
be defined as a function of t': -Empt^ 7 ) = vPbht'. Next, 
consider offloading in the interval (t',T], Fixed-rate trans¬ 
mission over this interval is assumed since this is the most 
energy-efficient data transmission policy under a deadline 
constraint as proved in |30|. Then the energy consumption for 
offloading, denoted as E 0 g(t'), can be written as E 0 ff(t') = 

L 2 

[2 — 1 ]2-(T — t'). The energy savings is thus given as 

[-UmptII 7 ) — E 0 ff(t')\. Varying t' changes the harvested energy, 
offloading energy consumption as well as energy savings. 
Specifically, as t' increases, Empt^ 7 ) grows linearly with t' but 
E 0 «(t') monotonically increases as shown in |30]. Therefore, 


the energy savings may not be a monotone function and 
thus optimization is necessary. To simplify notation, define 
the offloading duration t = T — t 1 and then EmfvU') and 
E 0 &(t') can be rewritten as Empt(A) = vPbh(T — t) and 
E a ff(t) = (2 at —1 )^-f. Substituting the expression of E MPT (f) 
and E a ff(t) rewrites the objective function as 

( 2 \ 2 
^-vPbhj t- <T ^t2Th. (14) 


Given this objective function, the problem for the current 
case is formulated as 

max E MPT (t) - E oS (t) 

s.t. 0 < t < T, (P4) 

EmptW — E 0 g(t) > 0. 
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2) Solution: Define 


p{h) 



In 2 

B 

i + wr-^f-\) 


(15) 


where W(x) is the Lambert function defined as the solution 
for W(x)e w ^ x ' > = x. Problem P4 is a convex problem as 
shown in the following lemma. 


Lemma 4 (Convexity of P4). The objective function of Prob¬ 
lem P4 is a concave function for t £ (0, oo) and maximized 
at t = p(h)L with p(h) defined in ( fl5] >. 

Proof: See Appendix [F] □ 

Define a positive constant a" as 


a 


n 




L In 2 
BT 




L In 2 
BT 


x exp + W{—e 1 J BT 2 ) + l^j 


(16) 


Then optimizing the objective function of Problem P4 over 
the interval (0, T) and investigating the feasibility condition 
yield the solution as shown in the following theorem. 

Theorem 2. The optimal offloading duration t* that solves 
Problem P4 satisfies the following. 

1 ) If Pbh 2 < a", the problem is infeasible. 

2) If Pbh 2 > a ", t* = p(h)L with p{h) defined in ( |T5| . 

Proof: See Appendix |G] □ 

Remark 3 (Maximum energy savings). It follows from The¬ 
orem [2] that the maximum mobile energy savings is 


^off = vPbhT - y[h)L 


where the function y{h) is defined as 


2/0) = a fjp 2 ex P ( W 


f vPbh 2 
V °’ 2 e 



(17) 


Remark 4. The tradeoff between the deadline T and BS 
transmission power Pb as discussed in Remark |T| for the 
local computing also holds for the current operation mode. 
Moreover, increasing the channel bandwidth B allows a more 
stringent deadline or smaller Pb- 

Remark 5 (BS transmission power vs. offloading duration). 
It can be observed from the expression of t* that increasing 
the BS transmission power Pb reduces the optimal offloading 
duration t*. The reason is that higher transmission power is 
affordable leading to a shorter transmission duration given 
fixed data to be offloaded. 


Remark 6 (Power beacon based MPT). A power beacon can 
be deployed for performing MPT such that the power transfer 
and offloading will be served by the power beacon and BS, 
respectively. This will lead to different gains for the MPT and 
offloading channels, instead of being identical in the current 
model, which, however, will not cause significant changes to 
the key results and policy structures. 


C. Offload or Not? 

Since the mobile has non-causal CSI, for each channel 
realization, it can decide the operation mode based on whether 
the successful computing conditions are satisfied and which 
mode achieves the larger mobile energy savings. 

First, if only one operation mode is feasible, i.e., h > a/Pb 
for local computing or h> \Ja"/Pb for offloading, then this 
mode is preferred. 

Next, if both operation modes are feasible, the desirable 
mode is selected by comparing the amounts of their maximum 
energy savings. Define the difference between their maximum 
energy savings as AS = S* ff —S* oc . It follows from Corollary [l] 
and Remark [3 that 

AS=^-y(h)L (18) 


where the coefficient 



and y(h) is given in (fTTb. Then offloading should be performed 

if and only if AS >V 

From m and m the effects of parameters such as the 
computation deadline T and BS transmission power Pb on the 
offloading decision are characterized as follows. 


1) If T < 


7 S 


y ( h j L , offloading is preferred which implies 
that a more strict deadline requirement tends to select the 
offloading mode. It can be interpreted as follows. As the 
deadline increases, the growing rate of harvested energy 
for local computing is larger than that for offloading. 
Moreover, the energy consumption for local computing 
decreases with an increasing deadline, however, that of 
offloading can be proved to be invariant from 


2) If Pb < -^(l + ea"' In a'") where 


7 0 


p/i 2 l 1TGU LLL “ ) VV11 ^ u eT 2 a 2 L In 2 ’ 

offloading is selected which indicates that offloading is 
preferred when the BS transmission power is insufficient. 
The reason for this threshold-based mode selection w.r.t 
the variation of BS transmission power is similar to that 
of deadline. 


IV. Energy Efficient Mobile Cloud Computing 
with a Dynamic Channel 

While wireless channel is assumed fixed in the preceding 
section, dynamic channel is considered in this section which 
is modeled as M channel fading blocks with block duration 
T c satisfying MT C = T. The mobile is assumed to have the 
prior knowledge of the channel power gains in these fading 
blocks (see Assumption [2}, enabling off-line data allocation 
for computing in different fading blocks. Similar to the fixed- 
channel counterpart, the mobile is assumed to select one of the 
two operation modes for a single task in these fading blocks. 
Local computing and offloading are optimized separately in 
the following sub-sections. The structures of resultant control 
policies are also analyzed. The results can be straightforwardly 
combined to optimize the mobile operation mode selection as 


in Section III-C with the details omitted for simplicity. 

















A. Energy Efficient Local Computing with a Dynamic Channel 

1) Problem Formulation: Exploiting the different time 
scales for channel variation and computing in Assumption [3] 
the problem of optimizing the CPU cycle frequencies for 
local computing can be decomposed into a master and a slave 
problems as follows. 

• Slave problem: Given a particular fading block and 
allocated data, the slave problem aims at minimizing the 
average energy consumption for computing in this block 
by controlling the CPU cycle frequencies in a similar 
way as the fixed-channel counterpart. Consider a single 
fading block with channel power gain h and allocated 
data size £. Moreover, let II > 0 denote the amount of 
residual energy from computing in the preceding block. 
The slave problem for this block is formulated as 


min 

{fk} 


s.t. 


M 0 

k =1 


'Yhlfl < R + vP b h^2 —, m =!,••• ,tN 0 , 


k=l 

£N 0 


k=l 


— u 


fc =1 


fk> o, k = l,---,£N 0 . 

(PS) 

Setting R = 0 reduces Problem P5 to PI for the case 
of static channel. Let Gi oc {£, R, h) denote the mini¬ 
mum average energy consumption for computing Gbit 
data in a single fading block. Then Gi oc (£, R, h) = 
7 Pfc(/fe) 2 where {/£} solve the above slave prob¬ 
lem. 

• Master problem: The master problem divides the input 
data for computing in different fading blocks under the 
criterion of minimum total energy consumption. Let n 
denote the index of fading blocks. The master problem 
is formulated as follows. 

M 

min ^ ] G\o C (£ n , Rni 6 n ) 

n= 1 

S.t. G\oc i Pn i hri) — Pn PvPbhnPci — 1? ? 1 


Pn — P-n—1 H"~ vPbhn—lTc 


G\oc^n— 1 ? Pn—lf hn— 1)5 


Pi — 0, 

£n> 0, 


M 

n—1 


n = 2, 


n = 1, 





Note that the first two constraints of Problem P6 imply that 
R n > 0 for n > 2. 

2) CPU-cycle Control Policy: The policy resulting from 
solving Problem P5 can be modified from that obtained from 
solving Problem PI for the case of static channel. Lor ease of 
notation, define two constants: 


b = 


f vP b hT g + RTl 

V 7^o 


n 


fvP b hT* + RTl 


7^1 


(19) 


and two energy consumption functions 

7V 5 o^ 3 


E 0 (£) = 


T? 


Exit) = 


TPiP 

T? 


where £ is the input-data size and do, 8i,(po and ipi are 
the scaling factors of £? for the two thresholds and energy 
consumption functions determined by the system model and 
CCI (see Appendix C of |18|), with 9q > 9\ and po < p\. 
Lollowing the same procedure as for deriving Theorem |T] and 
Corollary [l] the optimal policy for the current case is obtained 
as follows. 


Corollary 2. Consider an arbitrary channel fading block with 
the corresponding input-data size £ and residual energy R. 
The optimal CPU-cycle frequencies {fl } and the minimum 
average energy consumption G\ oc (£, R. h ) are described as 
follows. 


1) If £ < b , 


/ 1 ^ i\ _i 

fk = ^ p ™ ) p k 3 > and Gy.it, R - h ) = E o(t)- 

\ c m —1 / 


2) If b < £ < b'. 


ft = 


M 0 


iPm + A) 3 


m =1 


(Pk + A) 3 , Vfc 


where the positive constant A satisfies 


M 0 




Lfe=i 

and 


M 0 


+ A) 


.k -1 


vP b hT? + RTl 


Eq{£) < G\ oc (£,R,h) < E 1 (£). 


Moreover, when £ — > b, G \ oc {£, R, h) — > E 0 (£) and when 
t —> b' , G loc {£, R, h) — >■ Ei(£). 

3) If £ > b\ {/*,/a, • • • > fm 0 } is an em P l y set 


This Corollary shows that the slave problem is feasible only 
if £ < b' where b' is determined by the channel power gain h 
and the residual energy R. Moreover, when £ < b, the CPU- 
cycle frequencies and energy consumption are independent of 
Ph/i, implying that the energy harvesting constraint is inactive 
for small data-input size. However, for the case of b < £ < //, 
both {/£} and G\ oc {£, R,h) are determined by Pbh, which is 
consistent with the case of static channel. 

3) Sub-optimal Data Allocation Policy: The derivation for 
the optimal data allocation policy is intractable due to the lack 
of a closed-form expression for the energy consumption func¬ 
tion, Gtodi, R- h), which can be observed from the solution for 
the slave problem. In this sub-section, a sub-optimal but simple 
policy is derived by the approximation for Gi oc (£,R, h) and 
the amounts of residual energy in all fading blocks, denoted 
as {P„}. 

Lirst, the proposed sub-optimal data allocation policy re¬ 
quires only the following properties of G\ oc {£, R,h). 


Assumption 5. Lor an arbitrary block fading channel with 
data-input size £, Gi oc (£, R,h) is a monotone-increasing, dif¬ 
ferentiable and convex function for £ £ [0,6']. 
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The assumption can be justified as follows. First, the mono¬ 
tonicity of G\ oc (£, R, h) arises from the fact that computing 
more data requires additional energy. Next, finite computing 
energy per bit gives the rationality for the differentiability 
property. Last, the second derivative of Gj oc (£, R. h ) relates to 
the increase in energy consumption per additional data bit. For 
the special case of equal CPU-cycle frequencies, the energy 
consumption of T 2 per bit grows with data-input size t 9 
supporting the assumption on the convexity of G\ oc (£, R, h). 

Next, the residual energy variables {R n } are approximated. 
To this end, R n can be bounded as follows. 

Lemma 5. Given that local computing of input data is 
feasible, then R\ = 0 and 


ip(vP b h n -iT c + Rn-l) < Rn < vP b h n -iT c + Rn-1 (20) 


for n = 2, ■ ■ ■ ,M where tp = 1 — ^ > 0 . 

Proof: See Appendix |H| □ 

Note that the lower bound on R n is nonzero due to the 
energy harvesting constraint. Since it is difficult to obtain R n 
in closed form for the same reason as for deriving G\ oc (£, R, h) 
and the upper bound corresponds to the case without data 
computing, R n is replaced by its lower bound in the design of 
the proposed sub-optimal data allocation policy and the result 
is represented by R n . In other words, R„ = Cp(vP b h n _fT c + 
R n -i) for n> 2 and R\ = 0. Correspondingly, the constants 
b and b' defined in ( fl9| are modified as 

f _ ( v P h h n T'c + R n T’ r 2 \ 3 u _ ( vP b h n T*+R n T2\ 3 

"V 74 J “ " \ 7 ^i ) 


Using the above approximation and definitions, the minimum 
average energy consumption of the n-th fading block, denoted 
as Gioc, follows from Corollary [ 2 ] and Assumption [ 5 ] as 


G\oc[£ni Priori) — 



if t n < b n 
if b n < 4 < b' 


( 21 ) 


where <?(4) is a general function such that G\ oc (£ n , R n - h n ) 
has the properties in Assumption [5] 

Based on the above approximations, the data-allocation 
problem for minimizing the total energy consumption can be 
readily reformulated in a simple form as follows. 


(Sub-optimal Data Allocation) 


min 

t<r>} 

s.t. 


M 


Gloc(4> Rn 5 h n ) 


M 


n—1 

0 < 4 < v n , 


n = 1,2, ••• , M. 


(P7) 


Problem P7 is a convex optimization problem. The struc¬ 
ture of the resultant data allocation policy can be charac¬ 
terized as follows. Let &„(£) denote the root of equation: 

(b n , R n , h n ) = £ where £ is a Lagrange multiplier. Then 
the main result of this sub-section is obtained as shown below. 


Proposition 1. If L < Kt* Problem P7 is feasible. And 

the data-allocation policy that solves Problem P7 is given as 


l 


* 

n 


MO. 


h n 

h n > 


vP b T -3 

l<hbl(H)~R n T; 

vP b T | 


( 22 ) 


Proof: See Appendix [I] □ 

Note that {f*} are nonzero, indicating that the policy utilizes 
all fading blocks for computing since local computing and 
MPT can be performed simultaneously over all fading blocks. 


B. Energy Effcient Offloading with a Dynamic Channel 

1) Problem Formulation: Following the local computing 
counterpart, the problem for optimal computation offloading 
is formulated using the master-and-slave model as follows. 

• Slave problem: Given fixed allocated data, residual en¬ 
ergy and channel power gain, the slave problem aims 
at finding the optimal time division of a fading block 
for separate energy harvesting and offloading. Consider a 
single fading block with channel power gain h, allocated 
data-input size i and the residual energy R which comes 
from offloading in the preceding block. The slave prob¬ 
lem is formulated as follows for maximizing the energy 
savings in this block with the optimal time division. 

max E MVT (t, h) - E oS (t , h) 
s.t. 0 < t < T c , (P8) 

R + h) — E oS (t, h) > 0 

where S M pt(C h) = vP b h(T c — t ) and E 0 g(t, h) = (2b* — 
1 ) q jyt- Note that setting R = 0 reduces Problem P 8 to P4. 
Let G 0 a(£. R, h ) denote the maximum energy savings for 
offloading the Gbit data in this single fading block. Then 
G 0 fr(f, R , h) = Eyipj{t*, h) — E 0 ff(t*,h) where t* solves 
the above slave problem. 

• Master problem: The master problem concerns the 
optimal data allocation for offloading in different fading 
blocks with the objective of maximizing the total energy 
savings. Let n denote the index of fading blocks. Given 
the solution for the slave problem, the master problem of 
the optimal data allocation over different fading blocks is 
formulated as follows. 

M 

^ ' G 0 ff (4; Rn ■ bn) 

n—1 

n—1 

R n — ^ ^ ^off(^rri5 Rmi ^m )i Tl — 2, • , M , 

m—1 

R n > 0, n = 2, • • • , M, 

M 

Ri = o, Y, 4 = 4 

n—1 

4 > 0, n = 1,2, ••• ,M. 

(P9) 


max 

{C} 

s.t. 
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2 ) Optimal Time Division Policy: The slave problem. Prob¬ 
lem P8, can be modified from the fixed-channel counterpart. 
Problem P5, by adding the residual energy and thus solved 
following a similar procedure. For this purpose, define the 
following constants: 


T r B 


1 + W( 


vP b h 2 

a 2 e 


In 2 


and c = BT C log 



Then the optimal time division policy for the current case 
with residual energy R is obtained as shown in Corollary [3] 
following a similar procedure as for deriving Theorem [2] 

Corollary 3. Consider an arbitrary channel fading block with 
/:-bit data input and residual energy R. The optimal offloading 
duration t* and the maximum energy savings G 0 &{1, R, h) in 
this block are given as follows. 

1) If either: (a) R < - ^T c and i < vPb ^ +R or 

(b) R > BT {n 2 h ^ — a h T c and £ < c, then 


t* = p(h)£ and G oB (£, R, h) = vP b hT c — y{h)t 


where the constants y(h) and p{h) are defined in ( fl7j ) 
and respectively. 


2) If R > BT {^2^ — ar, d c 5= 5= c 'i then 

l <J 2 

t* = T c and G oB {£, R, h) = -(2^ - 1 )—T c . 

3) For other combinations of conditions for R and Prob¬ 
lem P8 is infeasible. 


In the above Corollary, the largest data-input size, optimal 
time division and maximum energy savings are determined 
by the channel power gain h and the residual energy II. 
Specifically, the time division of case 1) is the same as 
that of static channel case and the maximum energy savings 
decreases linearly with the growing of data-input size. Case 
2) corresponds to the scenario where the mobile has large 
residual energy and data input size. In this case, spending all 
time on offloading is the optimal time division policy. 

3) Sub-optimal Data Allocation Policy: One can observe 
from problem P9 that different summation terms in the objec¬ 
tive function are coupled due to the residual energy delivered 
from one block to the next. The conventional approach for 
solving this type of optimization problem is using dynamic 
programming (DP). The state space for the resultant DP is 
continuous and has to be discretized to facilitate iterative 
computation for the optimal policy, bringing high complexity. 
However, simulation reveals the sub-optimal low-complexity 
policy to be presented shortly can achieve close performance 
as DP. More importantly, the DP approach yields little insight 
to the structure of optimal policy while the said sub-optimal 
policy allows data allocation to be derived in closed form. 

The proposed sub-optimal policy is obtained by setting the 
residual energy variables as zero: {R n } = {0}, which is 
observed from the energy harvesting constraint in Problem 
P9 to be their lower bounds. Combining the approximation of 
{ R n } with Corollary [3] reduces Problem P9 as 


(Sub-optimal Data Allocation) 


max 

{«»} 

s.t. 


M 

^ ' vPbh n T c y[hn)£n 

n—1 


M 

n—1 


0 < In < 


vP b h n T c 

y{hn) 


(P10) 


n = 1,2, • • • , M. 


Problem P10 is a convex optimization problem and solving it 
gives the optimal policy in closed form. To state the policy, 
let the channel power gains {h n } be rearranged and re¬ 
denoted as {h n } such that {y{h n )} are in ascending order: 
y{hi) < y{h 2 ) < • • • < y(h.M)- Moreover, let n represent the 
permutation matrix with 


[hi,h, 2 , • • • , Hm] T = n x [hi, h. 2 , • • • , h M f. 


Using these definitions, the optimal policy from solving 
Problem P10 is given in the following proposition. 

Proposition 2. If L < ^f =1 vB \h" f c . Problem P10 is 
feasible. And it can be observed that the data-allocation policy 
solving Problem P10 is: 

\p* 0* . . . (/* ]T _ TT-1 v \p* p* ... p* ]T 

PmI 11 x 1*1 >-<-2) UmJ 

with {<?*} given in the following 

vPbhnTg 

ll(hn) ’ 

t _ sr^j vP b h k T c 

£-lk=l y(h k ) 

0, 

where j is unique and satisfies: 

vP b h n T c < ^ < X'' v Pbh n T c 

hi vCk) - hi yChn) 

The above data allocation policy is a greedy approach which 
allocates data to the fading blocks sequently by the ascending 
order of y(h n ) until all the input data has been allocated. 



n=l,2,---j 
n = j + 1, (23) 

otherwise 


V. Simulation Results 


In this section, the performance of wirelessly powered 
mobile cloud computing with static and dynamic channels is 
investigated by simulation. The parameters are set as follows 
unless specific otherwise. The data input size L is 1000-bit 
and the number of CPU cycles required for per bit is modeled 
by a Gamma distribution with a — 4 and /3 = 200 as in )T8| , 
resulting in e = 0.05. The constant determined by the switch 
capacitance 7 is 10~ 28 . The energy conversion efficiency v is 
0.8 GD- Let N t denote the number of BS antennas and set as 
N t = 2. The N t x 1 vector channel, denoted as h, is assumed 
to follow Rician fading and thus modeled as 


h = 


<IK - 

: h 


1 + K 


n 


1 + K 


where the Rician factor I\ £ {0,10}, the average fading 
power gain U = 5 x 10~ 6 , line-of-sight (LoS) component 
h has all elements equal to one, and h w is a Nt x 1 
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(a) Rician factor K = 0 


(a) Rician factor K = 0 


-f? 

f (4) /• 


* * 

* 

* 

iy 1 / 

/ / 1 5) 

— 

/ 


/ // 


/ / 1 

(1) MMS 

IJ' 

(2) Optimal local computing 

(3) Local computing with equal cycles 

(4) Optimal offloading 

(5) Offloading with equal partitions 


ol—- ' -'- 

0 0.05 0.1 0.15 

Deadline (s) 



(b) Rician factor K = 10 


(b) Rician factor K = 10 


Figure 2: Effect of deadline on the computing probability for the case 
of static channel. The BS transmission power is fixed at Pb = 0.5 W. 

i.i.d. CN{ 0,1) vector representing small-scale fading. The 
effective channel power gain h = ||h|| 2 , resulting from the 
transmit/receive beamforming. In addition, the variance of 
complex white Gaussian channel noise a 2 is 10 -9 W and the 
channel bandwidth B is 1 MHz. 

A. Static Channel 

The performance of three polices is evaluated, including 
the optimal local computing, optimal offloading and optimal 
mobile mode selection (MMS) integrating the preceding two 
polices. The baseline schemes for optimal local computing and 
optimal offloading are local computing with equal CPU-cycle 
frequencies and offloading with equal time partition (for MPT 
and offloading), respectively. 

Fig. [2] (a) and (b) show the curves of computing probability 
versus deadline T for the Rician factor K = 0 and 10, re¬ 
spectively. Several observations can be made. First, computing 
probability is observed to be a monotone-increasing function 
of T. Next, for a highly random channel (K = 0), the crossing 
of the curves (2) and (4) suggests mode switching as the dead¬ 
line increases. Specifically, local computing and offloading are 
preferred for relatively loose and strict deadlines, respectively. 
The reason is that compared with offloading, the computing 
probability for local computing grows faster as the deadline 
increases and also decays faster as the deadline decreases. Note 


Figure 3: Effect of BS transmission power on the computing proba¬ 
bility for the static channel. The deadline is fixed at T= 0.035 s. 



Figure 4: Approximation Gioc with = 1 W and T = 0.035s. 


that the thresholds for mode switching w.r.t the computing 
probability has no simple closed form. However, given LoS 
(I\ = 10), the optimal offloading is always preferred since 
the required transmission energy is small for such a channel. 
Last, compared with their corresponding baseline schemes, 
optimizing offloading shows more substantial performance 
gain than optimizing local computing. 

The curves of computing probability versus BS transmission 
power /), are plotted in Fig. [3] (a) for K = 0 and in 
Fig. [3] (b) for K = 10 with a fixed computing deadline 
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(a) Rician factor K = 0 


(a) Rician factor K = 0 



(b) Rician factor K = 10 


Figure 5: Effect of deadline on the computing probability for the 
dynamic channel. The BS transmission power is fixed at Pi, = 0.5 W. 



(b) Rician factor K = 10 


Figure 6: Effect of BS transmission power on the computing prob¬ 
ability for the dynamic channel. The computing deadline is fixed at 
T = 0.035 s. 


T = 0.035s. As observed from the figures, for a highly random 
channel (K = 0), the optimal local computing is preferred to 
offloading. However, for a channel with LoS (K = 10), the 
optimal local computing is preferred only when Pi, is small 
while optimal offloading should be chosen for large /',. Other 
observations are similar to those from Fig. [2] 

Consider the case of dynamic channel and data allocation 
over 4 fading blocks. The relatively small number of fading 
blocks is assumed to account for the difficulties of large-range 
channel prediction in practice. The sub-optimal data allocation 
policies derived in the preceding section are compared with 
two baseline polices: the optimal one computed based on the 
derived sub-optimal data allocation policy for local computing 
and DP policy for offloading, referred as DP data allocation, 
as well as a simple policy based on equal data allocation. 
Specifically, the DP policy for offloading is obtained by 
discretizing the state space and solving the Bellman equation 
backward recursively |[32). Each policy is feasible if one of its 
two operation modes is possible. For local computing, the min¬ 
imum average energy consumption for the n-th fading block 
Gi oc {£n, Rn, h n ) has no closed form. To allow simulation, the 
function with the properties in Assumption [5] is approximated 


as 


5(4) 


7<£(4)ln 

rp2 

C 


(24) 


where ip(£ n ) = (l - 4) 4 + <5o is a polynomial 

monotone-increasing and convex function satisfying <p(b n ) = 
ip 0 and — The approximation is verified by simula¬ 

tion to be accurate as shown in Fig. [4] 

Fig. [5] (a) and (b) show the curves of computing probability 
versus deadline T for the Rician factor K = 0 and 10, 
respectively. The proposed sub-optimal data allocation policy 
is found to have close-to-optimal performance and substantial 
performance gain over the equal allocation policy. Moreover, 
the gain for a highly random channel is larger than that for 
a FoS channel. This shows that adaptive data allocation is an 
effective way for coping with the effect of fading on mobile 
cloud computing. 

Fast, the curves of computing probability versus BS trans¬ 
mission power I), are plotted in Fig. [6] (a) for K = 0 
and in Fig. [6] (b) for K = 10 with a fixed computing 
deadline T = 0.035 s. Farge performance gain is observed for 
data allocation for the case of highly random channel. Other 
observations are similar to those from Fig. [5] 
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VI. CONCLUSION 

A novel framework of wirelessly powered mobile cloud 
computing has been proposed in this paper. Applying opti¬ 
mization theory, a set of policies have been derived for op¬ 
timizing the computing performance of two mobile operation 
modes, namely local computing and computation offloading, 
under the energy harvesting and deadline constraints. Further¬ 
more, given non-causal CSI, a sub-optimal policy for adaptive 
data allocation in time has been proposed to cope with the 
effect of fading on the system computing performance and 
shown to be close-to-optimal. This set of policies constitute 
a promising framework for realizing wirelessly powered and 
cloud-based mobile devices. 

This work can be extended to several interesting directions. 
First, full-duplex transmission can be implemented in the pro¬ 
posed system to support simultaneous MPT and computation 
offloading to improve the power transfer efficiency. Second, 
the current work focusing on a single-computing task can be 
generalized to the scenario of computing a multi-task program, 
which involves program partitioning and simultaneous local 
computing and offloading. Last, it is interesting to extend the 
current design for single-user mobile cloud computing system 
to the multiuser system that requires joint design of radio and 
computational resource allocation for mobile cloud computing. 


Appendix 

A. Proof of Lemma [7] 

Define the Lagrangian function for Problem P2 as 

N N / m m A 

^ = E i pk fk + E Xm E7/ 2 - vPbh E yk 


k =1 


m—1 


Kk=l 


k=1 


N 


N 


Vk 


(25) 


+ /W E» _T ) + E% ( 7T _: 

\fc=i / k =i 

Applying the Karush-Kuhn-Tucker (KKT) conditions gives: 


dL ( N \ 1 

^77 =27 Pkfk +27 fh f E Xm J-VkTT^ = 0, Vfc, (26) 

\m—k / 


dft 


U *) 2 


N 


M = - vPbh {]^f") + ^' lk=0 - 


A, 




Vk 




, fc =1 

N \ 

E^ T =0 ’ 

k =1 / 

1 


k= 1 


= 0 , 


ft 


-vt = 0, 


A fe > 0, p > 0, r) k > 0, fk > 0, y* k > 0, 

Then it is derived from ( |26| that 

p*\3 _ Vk 


Vk, (27) 
Vm, (28) 

(29) 

Vfc, (30) 
Vk. 


(.fir = 


27 (pk + ’Yhm=k 


Vk. 


(31) 


To ensure /£ > 0, since 7 and pk are positive and {A m } 
are nonnegative, it needs to satisfy that: rjk > 0 for all k. 


Combining it with ( [30] ) yields that yl = j*, leading to the 
desired result. ■ 


B. Proof of Lemma [2] 

First, due to the positivity of n k (see Lemma [TJ, it follows 
from ([27]) that: p. = rj k + vPbh Xk ) > Combining 

it with ([29]) gives: Ea^i Vk = Efc=i j* = T. 

Second, the optimal CPU-cycle frequencies © can be 
obtained by combining ( |27[ and ( |3 1 [ ). 

Third, compare /£ and f k+1 based on <[6]). Since p k +1 < p k 

and Em=fe+i < ELt A m , it follows that f* +1 > /*, 
completing the proof. ■ 


C. Proof of Lemma [j] 

It is proved by contradiction as follows. Assume there exists 
an integer m where 2 < m < N — 1 such that A m > 0. It 
follows from ( [28] ) that 


E 

k =1 
m— 1 

: E 


7 (flY - vP b h- 

Jk J 


k =1 L 


7 (/fe) 2 - vP b h- 
Jk J 


7(/m) 2 - vP bh — 

J m. 


= 0. 


Since ZZ=i 


7 (/fc) 2 — vP b hjr 

J k 


< 0, it can be obtained that 
7(/m) 2 ~ vPbh-pr > 0. Combining it with the monotonicity 
of {f k } (see Lemma |2] ), for the (m + l)-th cycle, it has 

7(/m+i) 2 - vP b h fi^ 1 > 7(/m) 2 - vP bhjv > 0, which 
results in 


m+1 r 


E 

k =1 
m 

E 

fc=i L 


7 (/fc) 2 - 

/fc 


7 (/ fc ) 2 - 

Jk J 


7(/m+i ) 2 -V p b h 


m+1. 


>0, 


contradicting the energy harvesting constraint. Similar proof 
by contradiction also applies to Ai. Therefore, the Lagrange 
multipliers {A m } for Problem P2 satisfy Ai = A2 = • • ■ = 
Ajv_i = 0 and A at > 0, yielding the desired result. ■ 


D. Proof of Theorem [7] 

Considering Problem P3, the conditions for feasible cases 
are derived as follows. 

1) Case 1: A > 0. First, substituting 0 into the deadline 
constraint leads to 


_V 


/i — vP^hX J 


N 


E+= + A )' 


,fc=i 


= T. (32) 


Combining (|7|i and ( [32] ) gives (|9). Next, substituting (|7]) 
into the energy harvesting constraint yields 


/ p — vP b hA \ 3 

V ^7 J 


" N 

E (Pk + A ) 


_k=1 


2 

3 


vP b hT 

7 


(33) 
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Combining ( |32| > and ( |33| results in ([TO]). Then, consider 
the monotone property of I),h. From ( [Tol l, the first 
derivative of P^h w.r.t A is 


d(P b h) 

dX 



where Uk 

w k = Pk 

that: 


t N \ 2 


- 


\k =1 / 

\k =1 

{Pk + A) 3, 

Vk 


VkW k 



X. Applying 


= (Pk + X) 3 ' 

Cauchy inequality 


(EfcLl VkWk 


■\N 


and 

gives 

. The 


(e£= iSt) > (Efe=i v k 

equality in (a) holds only when {^/ujfwk / y/vk/vJk} are 
equal for all k, which cannot be satisfied in this problem. 
Consequently, E 1 ’ < 0. Furthermore, the asymptotic 
properties of P b h, w.r.t A are characterized: when A —> 0, 
it has P b h — > a'\ when A — > 00, it has P b h —► a 
where a and a' are defined in <[8). Combining them with 
the monotone-decreasing property of P b h gives that: if 
0 < A < 00, it leads to a < Pbh < a'. In addition, when 
Pbh = a, there is only one feasible solution: /£ = ^ 
for all k, which can also be expressed as (|9]> by letting 
A = 00. 

2) Case 2: A = 0. Since is derived only using the 
deadline constraint, © can be obtained by letting A = 0 
in (|9]». Substituting it into the energy harvesting constraint 
(EfcLl 7 Uk) 2 ^ vP b hT) gives: P b h > a'. 

Thus, it can be concluded that if Pbh < a. Problem P3 is 
infeasible, completing the proof. ■ 


E. Proof of Corollary [7] 

The results of ( [T2| ), ( fl3] l and asymptotic properties can be 
derived straightforwardly following Theorem |T] The monotone 
property of E* oc for the case of a < P b h < a 1 is proved as 
follows. 

For notation simplicity, define Ofc = (pk + A)t such that 
Pk = af 3 - A. From © the first derivative of E^ oc w.r.t A is: 




By algebraic calculation, the part in the square bracket is 


N N 


£ £ E M ■- °5 + -«?«»)] 


i=l 

N N 


= £ £ E [<A(a? - a*?] > 0, 


i—1 j=i +1 


leading to > 0. Combing it with d( 'g b ^ < 0 and one-one 
mapping between Pbh and A results in that £j* c is a monotone- 
decreasing function of Pbh. ■ 


F. Proof of Lemma [7] 

2 

Define two constants for Problem P4: d - f — vP b h and 
d! = — . Then the first derivative of S 0 g w.r.t t for t G (0, 00 ) 

is given as 


3|s = «j + fd'-337l3') 2 #, 

dt \ Bt J 


(34) 


The second derivative follows: 


d 2 S oft d'L 2 (In 2) 2 l 

- = - —2 Bt < 0 

dt 2 B 2 t 3 


(35) 


since d! < 0, verifying the concavity of 5 0 ff. Then letting the 
hrst derivative © be zero gives: 


d 

dd 


L In 2 
Bt 



2 m . 


(36) 


Using the Lambert function, the solution for ( [36] l is: t = p(h)L 
with p(h) defined in ©• Furthermore, it can be observed 
from ( f34[ that if t —> 0, then —y 00 . Therefore, S a ff is 
maximized at t = p(h)L, leading to the desired result. ■ 


G. Proof of Theorem [2] 


Based on Lemma |4] if p(h)L > T, S* n is maximized at 
t = T and S'* ff < 0 such that Problem P4 is infeasible. 
Therefore, to guarantee the feasibility of Problem P4, two 
conditions should be satisfied: 1) t* = p(h)L < T; 2) 

s; s (t*) > 0. 

First, since —> 00 when t —> 0, it only needs to satisfy 

that when t = T, it has < 0. From ( |34| ), it can be obtained 
that 


r 


Pbh 2 > 


L In 2 
BT 


- 1 2 BT + 1 


(37) 


Next, substituting t* satisfying © into ( [T4| ) and letting 
S'off > 0 8 ives 


vP b h 2 L In 2 / / vP b h 2 

—exp [W ( — 
a z e B1 \ \ a z e 



(38) 


Denote d" = - \ and d!" = Then ([38]) is 

rewritten as d" + - > d'"e lv< - d "f 

e — 

Applying d" = W(d")e w( - d "' ) to the above in¬ 

equality and multiplying e~ d on both sides gives: 
\W(d") — d'"] e w ( d ")~ d '" > — e - l ~ d "\ Then, it follows that 


d" > 


W(—e 


— l—d 


0 d'" + W{-e~ 


(39) 


Substituting the expression of d" and d!" to ( [39] ) gives the 
solution for ( |38[ > as Pbh 2 > a" where a" is dehned in CD- 
Last, combine the two conditions ( |37| > and Pbh 2 > 

a". Since —1 < W(^~e~ l ~ d "') < 0, it has a" > 

In conclusion, P b h 2 > a" , com- 


(w-i) 2bt + i 


BT 

pleting the proof. 
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H. Proof of Lemma [5] 

According to Corollary [2] the mobile obtains the minimum 
average energy savings in this block if t n = b' n . In this case, 
the residual energy for the next fading block is 


Rn+ 1 — Rn + vPbh n T c — Gloe(C R ni h n ) (40) 

'yip, b 1 ^ 

= R n + vP b h n T c -P f f^. (41) 

Then substituting b' n given in ( fl9l ) into ( |4 1 [ i gives the lower 
bound. The upper bound is achieved when L n = 0. ■ 


I. Proof of Proposition [7] 

Define the Lagrangian function for Problem P7: 

M / M \ 

L = Y / G(en,Rn,h n )+dL-Y / Cn\ (42) 

n—1 \ n—1 ) 

M M 

T ^ ' W n ( in') + ^ ' C, n (in bfj • 

n= 1 n= 1 

Applying the KKT conditions leads to 


dL dG 
di n di n 
■w n i* n = 0, 


(C) Rn, h n ) - £ - W n + (n = 0, 

Cn (C - K) =0, > 0, Cn > 0, 


M 




Vra, (43) 
Vra, (44) 


First, it can be proved that i* n > 0 and w n = 0 for all n by 
the following steps: 

1) Observe from ( |2T| that r ?f {i* n , R n . h n ) > 0 and equals 
to 0 only when i* n = 0. 

2) Suppose there exists a n such that £* = 0. It leads to 
Cn = 0 and jff-(t n ,R n ,h n ) = 0. From <(43}, it gives 
C = -vj n < 0. 

3) There always exists one j where j f n such that i* > 
0. Therefore, Wj = 0 and (i*, Rj , hj) > 0. From 

( f43] i, it can be derived that Q = £ — (£*. Rj. hj) < 0 

which contradicts to the condition 0 > 0 and leads to 
the conclusion. 

Then, the data allocation follows: 

1) If t n < b’n > then Cn = 0 and §£(£*, Rn, K) = £• 

2) If £* = b’ n , then 

dC dC ^ 

— (t n ,Rn,hn) + tn= (i b' n , R n ,K) + Cn = £ 


In conclusion, i* n = min{6„(C), b' n } where b n (£) is the root 

of function §^-{b n , R n , h n ) = C and C satisfies J2n=i in = R- 

1 

Specifically, when b n (£) > b' n = ( vFbhn '^ c e + RnTc ^ 3 , it has 
h n < l6lbn ^pff ' nTc , completing the proof. ■ 
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